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1. Introduction 



One of the most important objectives of atomic 
physics in the last two decades, has been the achieve- 
ment of Bose-Einstein condensation in a weakly inter- 
acting gas. It is, therefore, not surprising that the obser- 
vation of Bose-Einstein condensation in the magnet- 
ically trapped alkali vapors "Rb [1] 'Li [2], and ^^Na 
[3], is considered to be one of the most exciting events 
of the past year. This is not in the least the case, because 
now that the successful road towards Bose-Einstein con- 
densation has been found, one can safely assume that 
gases of atomic *^Rb, '^^Cs, and 'H will also be Bose 
condensed in the near future and that a large number of 
different experimental systems will be available to study 
various aspects of the condensation phenomenon in 
detail. 

Moreover, both lithium and hydrogen have stable 
fermionic isotopes which can be trapped and cooled in 
a manner similar to that of their bosonic counterparts. 
For atomic ''Li this has already been achieved [4], but 
magnetically trapped deuterium (^H) has not been ob- 
served yet because the loading of the trap cannot be 
accomplished in the same way as in the case of atomic 



hydrogen. This is presumably caused by the fact that i) 
deuterium binds more strongly to a superfluid helium 
film, ii) the surface recombination rate is much larger, 
and iii) the sample is contaminated with atomic hydro- 
gen [5]. Nevertheless, it is to be expected that these 
experimental difficulties can be overcome and that one 
will soon be able to study both these Fermion gases in 
the degenerate regime. 

In view of the exciting developments mentioned 
above, it appears justified to present here an overview of 
the physical properties of all these weakly interacting 
quantum gases. We thereby intend to bring out most 
clearly the differences due to the Fermi or Bose statistics 
on the one hand, and due to the effective interatomic 
interaction being either repulsive or attractive on the 
other. This division in the focus of the paper is also 
reflected in its layout. In Sec. 2 we present a general 
introduction to the physics of the atomic quantum gases 
of interest. In particular, we specify the approximations 
that are allowed for these quantum systems, and subse- 
quently show how one can determine the effective inter- 
action between the atoms. Using the latter result. 
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we consider in Sec. 3 the Fermi and in Sec. 4 the Bose 
gases. Because we wish to bring out the different phys- 
ical properties most clearly, we treat there only the ho- 
mogeneous case. The experimentally more realistic case 
of a gas in a harmonic oscillator potential, however, 
does not lead to any qualitative differences, except for a 
Bose gas with an effectively attractive interaction. The 
latter possibility is therefore briefly discussed in Sec. 5, 
where we point out some problems that remain for the 
future. 

2. Interacting Quantum Gases 

All atoms mentioned in the introduction are members 
of the first primary group of the Periodic Table of 
Elements and, therefore, have a magnetic moment fi that 
is equal to the sum of the magnetic moments of the 
electron spin s (with s= 1/2) and the nuclear spin i . Due 
to this magnetic moment, the atoms can easily be 
trapped by an inhomogeneous magnetic field and can be 
evaporatively cooled by a resonant microwave field that 
flips the magnetic moment. The successful experiments 
of the past year have shown that this is a great advantage. 
However, the magnetic moment also gives rise to an 
important problem because two-body collisions in the 
gas can lead to a change in the component of the mag- 
netic moment parallel to the magnetic field and thus to 
a decay of the atomic density in the magnetic trap. 

This decay mechanism seriously limits the lifetime of 
the gas, and can in principle be caused both directly by 
the magnetic dipole-dipole interaction V^ and by the 
central (singlet/triplet) interaction V in combination 
with the hyperfine interaction V'^^ = a^fS ■ i [6] . In gen- 
eral, the typical time scale for the former process is 
much larger than that for the latter, and the gas sponta- 
neously spin-polarizes to a state in which the lifetime 
Tinei is determined solely by the magnetic dipole-dipole 
interaction. There are roughly speaking two scenarios: 
at small magnetic fields B « ahf/^te , where /tie denotes 
the electron magneton, the decay due to the central 
interaction can only be avoided by polarizing both the 
electron and the nuclear spin along the direction of the 
magnetic field, i.e., all the atoms are in the spin state 
\ms = 1/2, m, = ('). This corresponds to a so-called dou- 
bly spin-polarized atomic gas. However, for large mag- 
netic fields B » ahf/^te , the nuclear spin does not 
necessarily have to be polarized to obtain a relatively 
long lifetime and one can essentially trap any mixture of 
the states \ms = 1/2, m,) with an arbitrary projection of 
the nuclear spin. In this case we are dealing with an 
(electron) spin-polarized atomic gas. 

It is important to note that although the central inter- 
action does not determine the lifetime of the (doubly) 



spin-polarized gas, it does determine the typical 
timescale ta for elastic collisions. (See, however. Sec. 3 
for one exception.) Because ta « t^a for not too low 
temperatures, we can consider the spatial degrees of 
freedom of the gas to be in equilibrium, whereas the 
spin degrees of freedom clearly are not. For time scales 
short compared to the lifetime, we are thus allowed to 
apply equilibrium statistical mechanics to the gas and 
our main problem is the accurate treatment of the elastic 
part of the central interaction, which for the (doubly) 
spin-polarized gases of interest is the triplet interaction 

The determination of the effect of the interaction on 
the thermodynamic properties of a quantum gas is con- 
siderably simplified by the presence of two small 
parameters in the problem. In terms of the density n of 
the gas, the range of the interac tion ry , and the thermal 
de Broglie wavelength A = \/2Tr/i^lmk^T of an atom 
with mass m, the two small parameters are the gas 
parameter nry^ and the quantum parameter rylA. Physi- 
cally, the condition nry^ « 1 expresses the fact that we 
are dealing with a dilute system in which it is highly 
improbable that three particles are so close together that 
they can interact with each other simultaneously. As a 
result we can neglect three-body processes and only 
have to account for all two-body processes. Further- 
more, the condition r^lA « 1 shows that if two atoms 
interact, their relative angular momentum can only be 
equal to zero. In combination, we therefore conclude 
that for an accurate description of a quantum gas we 
only need to consider all possible two-body 5 -wave scat- 
tering processes. 

2.1 Two-Body Transition Matrix 

To explain how this may be achieved, we first of all 
introduce the appropriate "grand-canonical" Hamilto- 
nian of the (doubly) spin-polarized atomic gas. Using 
the language of second quantization [7], it reads 

H = 2j 2j i^k.a — l^a) ak,a ak,a 



"*" or/ 2^ 2^ ^q ^k+q,a tlk-q.o! ^k.a' tlk.a = Ho + Hmt, (1) 
aa' k'kg 

where ala creates and ak,a annihilates an atom in the 
(electron and nuclear) spin state I a) and the momentum 
eigenstate lA:) of a cubic box with volume y and periodic 
boundary conditions. As always, these creation and 
annihilation operators obey the (anti)commutation 
relations [a^.a, ai'>']± = 8k,k- S^.^ and [a^.a, a*>]± = 
[ala, at',a]± = if the atoms are fermions or bosons, 
respectively. Moreover, Vg = J dx Vjix) e.~"''' denotes 
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the Fourier transform of the triplet interaction and fi^ „ = 
fi^k^l2m + Sa = si, + Sa is the energy of the single particle 
states \k, a). Note that both the spin state la) and its 
energy s^ depend on the applied magnetic field B due to 
the Zeeman interaction V^ = - fi • B . For any particular 
atom this dependence can easily be determined by diag- 
onalization of the spin part of the atomic Hamiltonian, 
i.e., of y*^ + V\ but it is not explicitly calculated here 
since it plays only a minor role in the following. Note 
also that in agreement with our previous remarks, we 
have only assumed that the spatial degrees of freedom 
are in equilibrium and we have, therefore, introduced a 
different chemical potential jXa for each spin state la). 
Next we consider the effect of the interaction Hi„t on 
an "initial" state \i) of two atoms. In particular, we take 
\i)to be the state in which the occupation numbers of all 
the single particle states are zero except for Nxa+k.a and 
Nai-k.a; which are equal to one. In lowest order in the 
interaction, the energy of this state is Ei = sai+k.a + 
SK/2-k,a' if we use that the chemical potentials for this 
two-body problem are equal to zero. Applying standard 
perturbation theory, we find that this energy is shifted 
due to the interaction by an amount 



AEi- 






Vr-t + 






1 



Vk'-k+ 



A\k^-k'ym+iO 
which in the continuum limit y — > oo becomes 

V(k'-k) + . .. 



(5) 



1 



^Xk"- - k'ym+iO 



-8^,4 V{-2k) + J ^^ Y{-k-k') 



1 



^W-k'^)/m + iO 



V(k'-k) + 



(6) 



A£, = (il//i„JO + 



1 



2('-|//i„tl/);^— ^(/l//i„tlO + 



(2) 



f-^i 



where the "final" State I/) is an eigenstate of //o with 
energy Ef . Moreover, due to the possibility of scattering 
to other momentum states, the state \i) also acquires a 
finite lifetime T; , which up to second order perturbation 
theory is found from Fermi's Golden Rule: 



- = ^E 8(E-Ef)\{f\H;Ji)\- 



(3) 



Introducing the usual notation l/{x+iO) for the limit- 
ing procedure lim^io 1/ (x H- (17), both these results can 
be conveniently combined into 



The higher order contributions can be calculated in a 
similar way and we actually find that the complete per- 
turbation series can be summed byintroducing the so- 
called two-body T(ransition) matrix T^^{k',k; E), 
which obeys the famous Lippman-Schwinger equation 
[8] 



T'-\k',k;E) = V{k'-k) + 



/ 



dk" 



Vik'-k") 



1 



E-/i^k"^lm + iO 



T^\k",k;E) 



(7) 



and is proportional to the amplitude for two atoms to 
collide with a kinetic energy E in the center-of-mass 
frame and change their relative momentum from M to 
M'. Indeed, a comparison with Eq. (6) shows that the 
exact result is 



A£, - — = (i \H;Ji) + 2 ('■ l^int I/) 

2t, f 



{f\H,Ji) + 



Ei-Ef+ iO 



(4) 



The relevant matrix elements of the interaction Hi„t are 
easily evaluated and we obtain that 



AE-^ = ^\T'\k,k;A'kym) 



8a,^T^\-k,k;/i^kVm) 



(8) 



We therefore conclude that for a system of two atoms we 
can account for all two-body processes by first replacing 
in the Hamiltonian of Eq. (1) the potential Vq by 
the transition matrix T^{{k -k')/2 + q,(k -k ')/2; 
A^{k -k'f'/4m) and then using only first-order 



445 



Volume 101, Number 4, July-August 1996 

Journal of Research of the National Institute of Standards and Technology 



perturbation theory to avoid a double counting of the 
effects of the interaction. 

2.2 Many-Body Transition Matrix 

In a gas we are of course not dealing with only two 
atoms, but with a large number N=nV of atoms at a 
temperature T= Mk^^. This implies that we should re- 
peat our previous calculation for a grand canonical 
ensemble of "initial" states, in which the single particle 
states \k,a) have the average occupation numbers Nk,a = 
(exp;S(fiit,a - Ma) ± 1)' and the chemical potentials /Xa 
are determined from 



evident when we consider the contribution to the aver- 
age lifetime T; from a scattering process with a i= a'. 
Such a contribution is given by 



(1 -+- Nx/l+k-.a -+- NK/2~k;a) 

and consists of the sum of a loss rate 



2-17 1 V J^^k^ ^'^" 



N=^^Nk.a^^ N^. 



(9) 



I Vt-Vl' I Nxa +k,a Nki2 - k.d 



The average zeroth-order energy in this ensemble is 
equal to 



£■, = 2 2 if^k.a - Mo.) Nk,a 



(10) 



and its average energy shift and lifetime due to the 
interaction Hi^ can now be calculated perturbatively by 
averaging Eq. (4). Performing this thermal average is 
somewhat cumbersome, but ultimately we find 



i/i 



' aa Kk 



Nk/2 + k.a Niai-k.a' X | Vfl + — _2j 



,, 1 + Nk/2 + k\a + NkJ2 - k\a- ,, 



fi\k^-k'^)lm + iO 



(1 +NKI2^k;a) (1 +NKJ2~k;a) 

due to the collision of two particles, and a production 
rate 



277 1 V J^^k^ ^^k'^\ It. 12 

n V ^ \ m m 



Nxa+k.a Nk/2 ~ k.a' Nk/2+k;a A^B2-f,rf 

due to the collision of two holes. For a scattering process 
with a = a' we can obtain a similar interpretation, 
although it is slightly more complicated in this case due 
to the interference effects. 

Following our discussion of the two-atom problem in 
Sec. 2.1, we see that we can now include all two-body 
processes by introducing a many-body T(ransition) ma- 
trix r^f. (A:', k, K; E) which obeys 

TZ'{k\k,K;E) = V(k'-k) 



(K.....i2 



., 1 + NKa+k-,a + NK/2--k\a- ., 



(11) 



neglecting between the curly brackets terms that are 
either of order 0(^0 (nr^')'") or of order 0{VorvlA). 
(Diagrammatically speaking, we neglect the four 
"bubble" diagrams compared to the two "ladder" dia- 
grams.) 

The appearance of the factor 1 + NK/2+k\a- + NK/2--k\ct 
in this expression may appear surprising at first instance. 
It has, however, a clear physical meaning that becomes 



^f^^T/rfr' . „, 1+ JV.(g/2+k") + NAK/2-k") 
J (277)^ ^^^ " ' E-fi^k-^lm+iO 

TZ'{k",k,K-E) (12) 

in the continuum limit. Using Eq. (7), this is equivalent 
to 

TZ\k\ k, K; E) = T'-^k', k; E) 



.j^,T-ik;k',E) 



NaiK/2+k") + NAK/2-k") ^^s 



E-^^k'^-Zm + iO 



TZ-{k'\k,K-E), (13) 
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which is more convenient for our purposes since it offers 
the opportunity to explore the consequences of the 
small parameters nry^ and ry/A . Indeed, as a result of 
these small parameters the relevant momenta and ener- 
gies are always small compared to ^/ry and A^lmrv^, 
respectively, and we can safely neglect the momentum 
and energy dependence of the two-body T matrix. Con- 
sequently, Eq. (13) is immediately solved by 



T^. (k\k,K;E) = 

T^\0,0;0) 

1 ± r\0,0;0)Sa,aiK;E)' 



(14) 



where the quantity T^^{0, 0; 0) = ATraff'lm is propor- 
tional to the (triplet) i-wave scattering length a and 
Sa,a-{K; E) is defined by 



I 



i^a,a'(K; E) = 

dfc" NaiK/2+k") + NaiK/2-k") 



i27rf 



E-j^'-k"'-/m + iO 



(15) 



3. Fermi Gases 

We start our discussion of the interacting quantum 
gases with the Fermi case. We will only consider spin- 
polarized atomic gases in which two hyperfine states are 
almost equally populated. For convenience, these two 
hyperfine states will be denoted by IT) and \-l), respec- 
tively. The reason why we do not discuss doubly spin- 
polarized gases, is that in such gases the atoms can only 
scatter via ;? -waves due to the Pauli principle. As a 
result the (thermalizing) elastic collisions are dominated 
by the long-range magnetic dipole-dipole interaction V^ 
and the spatial degrees of freedom only equilibrate on a 
timescale that is comparable to the lifetime of the gas, 
i.e., Tei~Tinei- Furthcrmorc, if one is able to cool the 
system by other means, the gas basically behaves as an 
ideal Fermi gas on timescales short compared to the 
lifetime and interaction effects are negligible ([cf. Eq. 
(16)]). 

3.1 The Case a > 

We have argued that to account for all two-body 
5 -wave scattering processes we must use the interaction 



Moreover, using this solution we find that to a good 
approximation 



AEr 



2^ 



J_ 

2y 



2 2 NKI2+k.aNK. 



(1 + S,,,.) TZ{0,0,K; h^k'-lm). 



(16) 



which can again be reproduced by first-order perturba- 
tion theory if we replace V, in Eq. (1) by 

TZ\0,0, k + k'; /i\k - k'flAm). 

With the last conclusion we have basically achieved 
our objective of summing all two-body i-wave scatter- 
ing processes. In summary, we have shown that this 
requires the introduction of a many-body T matrix, 
which physically also incorporates the fact that the col- 
lisions of interest take place in a (gaseous) medium and 
are, therefore, influenced by the usual Fermi-blocking 
or Bose-enhancement factors [7]. As a result, the transi- 
tion matrix Ta,a'{k', k, K; E) acquires a dependence on 
the center-of-mass momentum ^K and the energy E that 
is explicitly displayed in Eqs. (14) and (15). This turns 
out to be important for the rest of the paper, where we 
use the many-body T matrix to arrive at an accurate 
(selfconsistent) theory of interacting quantum gases. 



^in. = ^Xt^,1 (0,0, k + k'; fi\k -k'flAm) 

' kk'q 

aLq.1 a^-q.l flf.J. a t,T. ( 1 7) 

However, this is not completely true because this inter- 
action induces on average also a change in the Hamilto- 
nian Ho, due to the fact that the occupation numbers 
Nk,a = {ak,a «t,a) do not vanish. This so-called mean-field 
correction is clearly given by 

AHo = 22 «*,««*:> j w 2 2 



Tff {0, 0,k + k'; h\k -k'flAm)Ni;,^ 



(18) 



and, denoting the density of atoms in the spin state la) 
by Ma, is well approximated by 



AHo — 2j 2j '^k.a O-k.a 
a k 

2 T'-\0,0-Q)nA, (19) 



since the dominant contributions to the integration 
over the momentum fik ' com e from the reg ion near the 
Fermi momentum M'p.a' = v 2m(^(,c.' - e^) (or from the 
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thermal momentum ^lA if the gas is not degenerate) 
where the many-body T matrix is almost equal to the 
two-body T matrix. The change in Ho can, therefore, be 
absorbed in a suitable redefinition of the chemical 
potentials. Introducing 



Ma = Ma- E T'^'{0,0-Q)n^, 



(20) 



mechanically stable. Indeed, using the hydrodynamic 
equations dnjdt = -V-Ja and Sa djjdt = - S/plm, we 
find that the two sound modes in the gas are described 
by the coupled wave equations 






m 



dpt/dnt dpiJdrii 
dpiJdm dpi/dni 



• V^ 



(24) 



the correct effective Hamiltonian for the spin-polarized 
Fermi gases finally becomes 



+ T^E T^'l(0,0,k+k';^\k-ky/4m) 

' kkq 



at+g,t<^t'-q,i^t',4-<^f,T ■ 



(21) 



where we have subtracted the constant VT^^(0,0;0) 
mni to avoid a double counting of the effects of the 
interaction that have already been included in the redefi- 
nition of the chemical potentials. Notice that we 
have implicitly assumed that the substitution 
jXa — > Ata is carried out also in the occupation numbers 
of Eq. (15), since this makes the calculation of the 
many-body T matrix self consistent. 

From this effective Hamiltonian we conclude that the 
spin-polarized gases are so-called Fermi liquids [9]. In 
particular, the elementary excitations (or quasiparticles) 
of the gas have a dispersion relation given by st^a - jJi^L 
and scatter off each other with an amplitude given by the 
many-body T matrix. As a result the density of atoms in 
the two spin states is determined by 



^a TT X J ^'fta 



(22) 



with the occupation numbers Nka equal to the Fermi 
distribution function (e'^^ + 1)' evaluated at £%,„ - /Xa. 
Moreover, the pressure of the gas is simply 



/'=¥EEln(l+e-^^--'^"') 

' a k 



+ mni = Zj Pa+ Pir 



(23) 



For a positive scattering length a , the interaction be- 
tween the atoms is effectively repulsive and the second 
term in the right-hand side is positive. Interestingly 
enough, this does not imply that the gas is always 



We thus conclude that the velocities of the sound modes 
are equal to the square root of the eigenvalues of the 
matrix in the right-hand-side of Eq. (24), and that both 
these eigenvalues need to be positive for the sound ve- 
locities to be real and the gas to be mechanically stable. 
In the case of a nondegenerate Fermi gas we can 
apply Boltzmann statistics to obtain pa = nJi-^T. The 
condition of mechanical stability then reduces to mnia^ 
^ {alA^IA and requires rather low densities. In the 
more interesting case of a degenerate Fermi gas, how- 
ever, we have 



riak^T 



Ha A' 



(25) 



neglecting terms of order 0{\l{kp^aA)'^). The condition 
on the densities then becomes ntnia'' ^ 'n-%304, which 
is much less restrictive. Note that in the mechanically 
unstable region, the gas is unstable against "spin- 
density" fluctuations and will phase-separate into two 
gaseous phases with the same total density mt + «j, but 
opposite "spin" densities nt - "J-- This process of phase 
separation is known as a spinodal decomposition and the 
density and temperature conditions at which it first oc- 
curs, lie in our case on a so-called spinodal surface. 
Apart from the first-order phase transition implied by 
the above instability, we do not expect any other phase 
transition to occur in the gas on the basis of our effective 
Hamiltonian. In particular, we do not expect a phase 
transition due to quantum degeneracy, because for that 
we need an attractive interaction between the quasiparti- 
cles, as will be explained shortly. In principle, however, 
it is possible that a positive i-wave scattering amplitude 
induces a negative p -wave scattering amplitude due to 
the Kohn-Luttinger effect [10]. In this manner a phase 
transition to a superfluid state might occur after all. 
Although it is interesting to analyze this possibility in 
more detail, we will not do so here since it is beyond the 
scope of the present paper. 

3.2 The Case a < 

The prospects for achieving a phase transition in the 
degenerate regime are more promising for spin-polar- 
ized Fermi gases with a negative scattering length, 
which actually applies to both atomic deuterium [11] 
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and atomic ^Li [12]. Nevertheless, for this phase transi- 
tion to be experimentally observable, it must take place 
in the (meta)stable region of the phase diagram where 
the gas is mechanically stable. Since the stability analy- 
sis of the previous section is insensitive to the sign of 
scattering length, this implies in the degenerate regime 
that the densities of the two spin states have to fulfill the 
requirement ntnia'' ^ TrV2304. In contrast to the case of 
a Fermi gas with a positive scattering length, the latter 
condition now does not rule out the possibility of the 
formation of a superfluid state at sufficiently low tem- 
peratures. 

The physical reason for this is that with effectively 
attractive interatomic interactions the gas can form a 
Bose condensate of Cooper pairs, in the same way as an 
electron gas can form a condensate of Cooper pairs in 
the BCS theory of superconductivity [9]. Mathemati- 
cally, the instability of the gas towards the creation of 
Cooper pairs is signalled by the fact that the scattering 
amplitude diverges for two quasiparticles at the Fermi 
energy Sp = /i^kll2m = {/i^kl^^l2m + h^ klJ2m)l2 with 
opposite momenta and "spin" [13]. Using Eqs. (14) and 
(15) and neglecting the small imaginary part of 3{0; 
2e-e), we thus find that the critical temperature of the 
BCS transition obeys 



1 



{0,0,0; 2 Sp) Airah- 



+ P 



Ak 



N,(k) + N,(k) _ Q 
2sp - fi^k^lm ' 



(26) 



where p denotes the Cauchy principal-value part of the 
integral. 

The properties of this so-called linearized BCS gap 
equation are well-known from the theory of supercon- 
ductivity. First, this equation only has a solution if a < 0. 
Second, the highest critical temperature for a given total 
density n is obtained if nt = «i or equivalently if 8ep = 
\/i^kl.^l2m -A^klJ2m\ = 0. This can be understood 
most easily by noting that for a nonzero value of 8ep 
there are less particles at the Fermi surface that can be 
paired, which suppresses the critical temperature. 
Third, if the densities in both spin states are equal, the 
gap equation can be solved analytically [14], and leads to 



"spin polarization" Int - «il/(«T + «i) is less than SkgTJ 
2sp [12]. Qualitatively, this is a result of the fact that if 
Sfip is of the order of ksT^, it is no longer energetically 
favorable to form Cooper pairs and the gas is in a normal 
state even at zero temperature. 

From an experimental point of view, we thus con- 
clude that the most favorable conditions for the observa- 
tion of the BCS transition in a spin-polarized Fermi gas 
are obtained if nt = ni = n/2. Under these conditions the 
mechanical stability of the gas at the temperatures of 
interest requires that n la P ^ Tr/24 or equivalently that 
kp\a I ^ tt/2. For atomic deuterium the triplet scattering 
length is equal to a — -6.8 Oo (flo is the Bohr radius). 
At realistic densities of ~ 10'^ cm~^ we thus find that 
kp\a\ — 1.1 X 10' ^ « tt/2 and the stability condition is 
easily fulfilled. However, from the exponential depen- 
dence in Eq. (27) we immediately see that the critical 
temperature is completely out of reach with present-day 
cooling techniques. Fortunately, this conclusion does 
not hold for atomic ^Li because the triplet scattering 
length now has the anomalously large value of 
a — -4.6 X 10^ flo- For the same density of atoms we 
thus find that kp\a\ — 0.4 < tt/2 and that the critical 
temperature is only 29 nK, which is rather close to the 
temperatures that recently have been obtained with the 
bosonic isotope ^Li [2]. Therefore, it appears that ''Li is 
a very promising candidate for the observation of the 
BCS transition in a weakly interacting gas. 

Because of this possible application, we also want to 
mention briefly how we can arrive at an accurate de- 
scription of the gas below the critical temperature. 
Although it is not difficult to treat the general situation 
with nr i= rii, we concentrate here on the most impor- 
tant case of an equal population of the two hyperfine 
states. Moreover, we also only discuss the derivation of 
the elementary excitations of the superfluid state (which 
is sufficient for the most interesting equilibrium proper- 
ties of the gas) and do not consider the effective interac- 
tion between these quasiparticles. Introducing the con- 
venient notation jx' = jXa - Sa for the renormalized 
chemical potential (or Fermi energy Sp) of both spin 
states, our task is therefore reduced to including into the 
free Hamiltonian of the quasi particles, i.e., into 



^" = E E (fit - At') «M«t,. -VT'^'iO, 0; 0) J, (28) 

a k 



T ~ ^ 

^' SkB^'^Pl 2kp\a\ 



(27) 



For unequal densities an analytical treatment is not fea- 
sible, but a numerical solution of Eq. (26) shows that 
nonzero critical temperatures are only possible if the 



the effect of the Cooper pairs. 

To achieve this, we must realize that BCS theory is 
roughly speaking the theory of a Bose-Einstein conden- 
sation of Cooper pairs. Therefore, the order parameter 
of the phase transition is the expectation value of the 
annihilation operator for a pair of atoms. (See Sec. 4 for 
a treatment of Bose-Einstein condensation in an atomic 
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Bose gas.) As a result, we do not only have nonvanishing 
average occupation numbers {ai^a «t,a ) but also nonvan- 
ishing values for {ak,a fl-t.a)- Note that in our case we are 
dealing with a pairing between atoms with different 
spin states and the latter average is only nonzero for 
ai= a'. To see what the influence of the anomalous 
averages {ak,fa.k,i) is, we briefly return to the Hamilto- 
nian of Eq. (1). Taking for the sake of clarity V, = Vo, 
the lowest order average correction on the Hamiltonian 
Ho due to the interaction now becomes [cf. Eq. (19)] 

a k ^a' ^ a ■' 

+ 2 {^oflt.T, aXi + Aoa.k.i at,?} , (29) 

k 

if we define the BCS gap parameter Ao by 



^ k 



(30) 



From our previous discussion of the normal state of the 
gas we know that we can to a good approximation 
account for all two-body processes by replacing in 
Eq. (29) the potential Vo by r"'(0 , ; 0). It is important, 
however, that we do not perform this substitution in Eq. 
(30) because the BCS theory is already going to sum all 
the other relevant two-body processes for us [15]. 

We thus arrive at the conclusion that the free Hamilto- 
nian of the quasiparticles, including the effects of a 
condensate of Cooper pairs, is given by 

^o"'' = 2 2 (s* - m') ala ak.a + 2 



[Aoalt alkA + Al a_t,j,at,T} 



y^-vr^«(0,0;0)^, 

Vn 4 



(31) 



where we have again subtracted a constant term to avoid 
a double counting of the interaction effects. Being 
quadratic in the creation and annihilation operators, 
this Hamiltonian can be diagonalized by means of a 
Bogoliubov transformation, and we find that the Bogoli- 
ubov q uasiparticles ha ve an energy dispersion 
^&)t = V (fit - fx'Y + l/ioP which for nt = mj- is indepen- 
dent of their "spin." Moreover, after this diagonalisation 
the average density n and order parameter Ao are easily 
calculated and result in the equation of state 



and the BCS gap equation 



— -Y 



l-Nk,t-Nk,i 



ih 



U>k 



(33) 



respectively. Here, the average occupation numbers of 
the quasiparticles are again denoted by Nk.a , but they are 
now equal to the Fermi distribution function evaluated at 
Auiic and are in fact independent of a. 

The equation of state and the BCS gap equation are, 
at a given density and temperature, two equations for the 
two unknown quantities /a' and \Ao\. However, before we 
can actually solve these equations we first have to resolve 
the ultraviolet divergence in the gap equation. This diver- 
gence is a consequence of the fact that we have 
neglected the momentum dependence of the potential 
and used V, = Vo . Making use of the same approxima- 
tion in the Lippman-Schwinger equation Eq. (7) and of 
the fact that the Fermi energy /a' = A^k^llm is always 
much smaller than fp'tmry, we find that the two-body T 
matrix obeys 



1 



1 



1 



Is 



1 



T'^^{0,0;0) T^''(0,0;2ijl') Vo V^2(st 



(34) 



in this case. We thus observe that the divergence in the 
BCS gap equation can be cancelled by a renormalization 
of l/Vo to l/r"'(0,0;0) = m/47ra^\ In this manner we 
then obtain the following gap equation 



4TTa^' 



+^2 

* k 



l-Nk,t-Nk,i 



1 



2(et - ix') 



(35) 



which is free of divergences. Furthermore, taking the 
limit /io — > to obtain an equation for the critical tem- 
perature Tc, we exactly reproduce (after taking also the 
thermodynamic limit and substituting /a' = Sp) the result 
of Eq. (26) which shows the consistency of our approach 
and, in particular, also that the BCS theory in essence 
indeed sums all two-body scattering processes. 

4. Bose Gases 

We now turn our attention to the atomic Bose gases, 
which have caused so much excitement in the atomic 
physics community after the experimental observation 
of Bose-Einstein condensation in atomic *'Rb, 'Li and 
^^Na vapours. As in these experiments, we consider here 
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only the doubly spin-polarized case. This implies that 
there is only one spin state in the problem and we will, 
therefore, in this section completely suppress the spin 
indices in our notation. 

4.1 The Case a > 

Again, we first consider a gas with an effectively 
repulsive triplet interaction, which applies for example 
to atomic hydrogen [6], atomic ^'Rb [16] and atomic 
^^Na [17]. The appropriate effective Hamiltonian for the 
doubly spin-polarized atomic Bose gases is, in analogy 
with the discussion in Sec. 3.1, given by 

//<=« = 2 (fit - M ')«*«* - VT'^O, 0;0W 



want to describe the gas below the critical temperature, 
we need to consider the influence of a nonvanishing 
average (ag) on the effective Hamiltonian. Focussing 
again on the derivation of the elementary excitations in 
the superf luid state, this is most easily achieved by first 
carrying out the substitution a^ -^ S^o (a^) + a^ in the 
Hamiltonian H"'', and then keeping only all terms that 
are of first or second order in the creation and annihila- 
tion operators. In this manner we find in a good approx- 
imation that 



Hf = -V r\0 ,0;0)(n 'f + 2T'^\0 ,0,0;0)n'no 



+ ^T''\0,0,0;Q)nl 



+ :r^2 T''\0,0,k + k';/i\k-k'flAm) 

^" kkq 



O-k + qO-k'-q O-k- O-k, 



(36) 



with a renormalized chemical potential jx' equal to 

IJL' = IJL-2nT^\0,0;0). (37) 

Looking at this hamiltonian we immediately see that the 
density of atoms is determined by 



-k|"' 



(38) 



where the average occupation numbers Nk are now equal 
to the Bose distribution function (e'^'' - 1)"' evaluated at 
fit - /a'. In addition, the pressure of the gas is 



IcbT 
V 



Eln(l 



-/3(q-rt-) + '^■^a^ ^2 



(39) 



which for a gas with a positive scattering length does not 
show a mechanical instability, i.e., dp/dn is always 
greater than zero for densities such that nry^ « 1 and 
our theory is applicable. 

On the basis of Eq. (38) we thus expect the gas to 
Bose condense at the same critical temperature 



To: 



2ttA^ 
mks V^(3/2) 



(40) 



as the ideal Bose gas. For an interacting Bose gas with 
positive scattering length, the order parameter corre- 
sponding to this phase transition is the expectation value 
of the annihilation operator for an atom [18]. If we also 



-lJL' + noT'"'{0,0,0;0)][ao + a^ 



+ 2 («t - M')«tV + E [(2noT''\0, 0,0; 0))at^at 
+ 1 Mo T^'^O, 0,0; 0) (aU\ + fl-kflt)}, (41) 



if we choose {ao ) to be real and introduce the condensate 
density no = {cioYlV and the density of noncondensed 
atoms n' = n - Mq- The most difficult part of the calcula- 
tion is again to avoid double countings of the interaction 
effects. Besides the subtractions in the first term 
in the right-hand-side, this also requires that the relation 
between fx and fx' is changed into iJi = 2n' 
T^\0 ,0;0) + ix', which reduces to Eq. (37) if the con- 
densate density vanishes. 

To proceed, we first note that if we want (ag ) to be the 
total expectation value of the original operator ag, we 
must require that the expectation value of the new oper- 
ator ao (i.e., after the substitution) disappears. This is 
achieved by putting ix' = no T'^^{0,0,0;0), because 
then the terms in the hamiltonian linear in the creation 
and annihilation operators vanish. As a result, the chem- 
ical potential is determined by 

ix = 2{n- no) T^%0 ,0;0) + no T'^^O ,0,0;0), (42) 

which for the approximation that we are using is just the 
famous Hugenholtz-Pines relation [19]. Because 
our theory obeys this relation, we expect that the disper- 
sion relation A wt of the quasiparticles is linear at long 
wavelengths. This indeed turns out to be the case, since 
a diagonalization of the above hamiltonian by means of 
a Bogoliubov transformatio n shows that ^cot 
= ^e\ + 2noT''\0,0,0;0)e, [20] and that 
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Hf=J,Aw,blh 



V\ r\0, 0; 0)(«') + 2T'"'iO, 0,0;0)n'm 



+ ^T'"'{0,0,0;OW], 



(43) 



where bl and bu denote the creation and annihilation 
operators for the Bogoliubov quasiparticles, and we have 
been neglecting between the brack ets additional terms 
that are a factor of order 0(yno(P) smaller [18]. More- 
over, using this Bogoliubov transformation we can now 
obtain also for the superfluid state the pressure 

+ T^\0, 0;0)(n'f + IT^^O, 0,0;0)n'no 
1 



+ ^T'"(0,0,0;0)no^ 



(44) 



and the equation of state 

M = «0 + 2j\ ^ Nk 



st + uo T'^\0,0,0;Q)- hci)t 



2^u)k 



(45) 



which determines the condensate density at a fixed den- 
sity and temperature. Note that in the equation of state 
the average occupation numbers of the quasiparticles Nk 
are equal to the Rose distribution function evaluated at 

This almost completes our discussion of the interact- 
ing Bose gas with a positive scattering length. However, 
there remains one point that we need to address, namely 
that our theory is not yet self consistent. This is because 
we have not yet mentioned how the many-body T matrix 
is determined below the critical temperature T^. In gen- 
eral, this is a complicated problem due to the presence 
of infrared divergences in the theory of a dilute Bose gas 
[21]. At very low temperatures such that naA^ » 1 the 
treatment of these divergences requires more advanced 
methods then the ones we are using here. Therefore, we 
consider only the opposite limit naA^ « 1, which is 
in fact the most relevant one for experiments. In this 
regime the average kinetic energy k^T of the atoms is 
much larger than the average interaction energy 
ATTan/P'lm and the Bogoliubov dispersion is for 



all practical purposes well approximated by 
/ia)^ ^-^ s^ + no r*"'(0, 0,0;0) = Si, + /i,'. In Eq. (41) this 
means that we can neglect the anomalous terms propor- 
tional to alaJ^ or a^ai^. Hence, we now basically deal 
with the Hamiltonian 2(e;i + /x') a^aic and the many- 
body Tmatrix is again just given by Eqs. (14) and (15). 
Finally, we have to deal with one last subtlety. In the way 
we have defined the many-body T matrix, it describes 
the scattering of quasiparticles. However, we are dealing 
with a Bose condensate of particles and we actually 
need their scattering amplitude. In our case, fortunately, 
this only leads to a slight shift E -^ E-lfi' in the en- 
ergy of the many-body T matrix. We thus arrive at a 
fully self-consistent theory of the superfluid phase in 
the regime naA^« 1 , if in all formulas below Eq. (40), 
T'^^(0 , 0,0; 0) is replaced by the many-body T matrix 
r^''(0,0,0;-2;u,'). 

The transition matrix r^^(0, 0, 0;-2/x') is almost 
equal to T^^(0, 0;0) = ATTZ./f'lm for most temperatures 
below Tc, but decreases to zero for temperatures near the 
critical temperature due to the fact that H(0;-2(x') 
diverges in the limit /x' — > 0. Therefore, our self-consis- 
tent approach reduces to the conventional Bogoliubov 
(or Popov) theory of a weakly interacting Bose gas for 
most temperatures. However, near the critical tempera- 
ture it resolves a long-standing problem of the Bogoli- 
ubov theory, which actually predicts a first-order phase 
transition [22] instead of the second-order phase transi- 
tion expected from the theory of critical phenomena. Of 
course, using the many-body T matrix instead of the 
two-body T matrix does not include all the critical fluc- 
tuations that are of importance near the critical temper- 
ature. Nevertheless, using the many-body T matrix we 
incorporate, at least qualitatively correctly, the fact 
(known from renormalization group theory [23]) that 
the interaction between the condensate particles should 
vanish exactly at the critical temperature. It is rather 
interesting that this can already be achieved by including 
only all two-body scattering processes [24]. 

4.2 The Case a < 

We now turn to the physics of a Bose gas with effec- 
tively attractive interactions, which is realized for 
instance in atomic 'Li [25] and atomic '^^Cs [26]. This 
is an interesting topic because a thoughtless application 
of the theory of Sec. 4. 1 would result in a dispersion of 
the Bogoliubov quasiparticles that is equal to 
hwk = Ve| - (8 77 1 a I Uoff' Im ) fii and therefore purely 
imaginary at long wavelengths. From this simple argu- 
ment we can thus already conclude that Bose-Einstein 
condensation in the canonical sense, i.e., with the order 
parameter {ao ), will not occur if the scattering length is 
negative [27]. However, from the experience with 
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fermionic gases gained in Sec. 3.2, we might expect a 
BCS-like phase transition to take place instead. This 
would be signalled in Eq. (36) by a divergence of the 
scattering amplitude of two quasiparticles at energy /a' 
and with opposite momenta, or by [cf. Eq. (26)] 



1 



T''\0,0,0;2iJ,') 



1 



r%0,0;0) 



f dk 

J (277 



Nik) 



y iJi' - fi^k^llm 



= 



(46) 



with of course N{k) the Bose distribution function eval- 
uated at fi^k^l2m - /a'. Analyzing this condition, we 
find that the phase transition indeed occurs at a temper- 
ature Tbcs = 7o(1 -h 0( |a |/Ao)), which is only slightly 
above the critical temperature of the ideal Bose gas. 

Following the same procedure as in Sec. 3.2, we can 
easily show that below this critical temperature the 
equation of state becomes 



-i?{ 






Nu + 



e/c - fJi' - ^cok 



(47) 



and that the BCS gap equation is well approximated by 



1 



Aiiafi^ V J ^(O/c 



1 



Nk 



(48) 



Here the energy of the Bogoliubov quasiparticles is 
^wt = V(fiit - fJb'y - \Ao P and their average occupation 
numbers Nt are equal to the Bose distribution function 
evaluated at this energy. These equations again deter- 
mine jx' and \Ao\ at a fixed density and temperature. 
Notice that in contrast with the fermion case, the disper- 
sion ^cok has a minus sign in front of Uo P- This has 
important consequences because it implies that if 
we lower the temperature at a fixed density, both Uo 1^ 
and jx' increase in such a manner that the gap in the 
Bogoliubov dispersion decreases. Hence, at a second 
temperature Tbec < Tbcs, the gap closes and the number 
of particles in the zero-momentum state diverges, which 
signals a Bose-Einstein condensation. Below that sec- 
ond critical temperature we always have \Ao I = - ^t ' and 
the dispersion ^co^ becomes equal to \/e^ -liju'e^ 
which, interestingly enough, is exactly the same Bogoli- 
ubov dispersion as for a Bose-condensed gas with posi- 
tive scattering length. Moreover, Eqs. (47) and (48) turn 
into 






Sk- IJ' 

hu)k 



Nk + 



Bk - fJi' - filSik 

2hu)k 



(49) 



Airah^ 



M 



V 



1 



/iwk 



(50) 



and 



determining now ^t' and the condensate density no- It is 
important to mention here that the possibility of a Bose- 
Einstein condensation in a gas with effectively attractive 
interactions, in the way that we have just seen, does not 
contradict the argument against Bose-Einstein conden- 
sation presented in the beginning of this section. The 
reason is that in that argument we assumed that a Bose- 
Einstein condensation would be associated with the or- 
der parameter {uo) and did not consider the possibility 
that it would be associated with the condensate density 
Mo, as in the ideal Bose gas [27]. 

It thus appears that a Bose gas with a negative scat- 
tering length offers the exciting possibility of observing 
both a BCS transition and a Bose-Einstein condensation. 
However, we have up to now not considered the mechan- 
ical stability of the gas in these two superfluid phases. 
This is clearly of the utmost importance, because Eq. 
(39) shows that for a negative scattering length, dp/dn 
does not always have to be positive. Indeed, calculating 
the spinodal line from dp/dn = 0, we actually find the 
unfortunate result that the gas is always unstable at Tbcs- 
We are thus forced to conclude that the two second- 
order phase transitions that we have found above are 
always preempted by a first-order gas-liquid or gas-solid 
transition and are, therefore, in a homogeneous system 
unobservable. 

5. Discussion and Outlook 

In this paper we have presented a unified picture of 
weakly interacting atomic gases, focussing on the differ- 
ences that arise due to the statistics of the atoms and due 
to the interatomic interaction being either effectively 
repulsive or attractive. However, we have restricted our 
discussion to the homogeneous case, whereas the exper- 
iments are up till now always performed in an inhomo- 
geneous magnetic trap. Fortunately, the experiments are 
also always performed at temperatures such that 
kBT» ^M, where ^w is the energy splitting of the (har- 
monic oscillator) trapping potential. Under these condi- 
tions we are generally justified in using a local-density 
approximation, which means that we consider the gas as 
homogeneous in each point in space. In this sense the 
theory presented above can also be applied to an inho- 
mogeneous situation and, in first instance, we do not 
anticipate any qualitative differences. Nevertheless, for 
quantitative predictions of, for example, the critical tem- 
perature and the density profile in the trap, it is of course 
essential to include the effect of the trapping potential. 

Moreover, for a Bose gas with negative scattering 
length there is also a qualitative difference. In Sec. 4.2 
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we argued that for a homogeneous gas a Bose conden- 
sate is always unstable. However, first Hulet [28] and 
subsequently Ruprecht et al. [29] have argued that in an 
inhomogeneous situation one can have a metastable con- 
densate if the number of particles is sufficiently small. 
This was recently put on a more firm theoretical basis 
by Kagan et al. [30] and by one of us [3 1]. In particular, 
the last author has shown by a calculation of the decay 
rate of the condensate due to thermal fluctuations, that 
a long-lived metastable condensate can exist for temper- 
atures sufficiently close to the critical temperature. Al- 
though this result appears to explain the recent observa- 
tion of Bose-Einstein condensation in atomic ^Li, our 
theoretical understanding of these experiments is still 
incomplete. The most important problem that remains to 
be solved is the stability of the gas just above the critical 
temperature. Clearly, on the basis of a local-density 
approximation we would conclude, as in Sec. 4.2, that 
the gas is unstable, and that also in an inhomogeneous 
situation Bose-Einstein condensation is unobservable. 
However, the local-density approximation is invalid near 
the critical temperature due to the large value of the 
(homogeneous) correlation length. A thorough discus- 
sion of this point, therefore, requires a theory that goes 
beyond this approximation and is evidently of great in- 
terest in view of the Rice experiments with ^Li [2] . 

Also in the case of a Bose gas with positive scattering 
length, interesting theoretical problems remain. One of 
them is the correct treatment of the infrared divergences 
in the regime naA^ » 1, which, due to the effectiveness 
of evaporative cooling techniques, is also anticipated to 
be experimentally accessible. As already mentioned, we 
presumably need more advanced methods than dis- 
cussed here to accurately perform a resummation of the 
divergent contributions to the thermodynamic quantities 
of interest. Indeed, we are presently carrying out a renor- 
malization group study of the interacting Bose gas in 
this regime, because by this method infrared diver- 
gences are automatically avoided. Moreover, a resolu- 
tion of these divergences is also of importance for the 
theory of two-dimensional atomic Bose gases, such as 
doubly spin-polarized atomic hydrogen adsorbed on su- 
perfluid "^He films. The reason for this is that in a two- 
dimensional geometry, the divergences are so important 
that they actually prevent the possibility of a Bose con- 
densate (in the thermodynamic limit). Nevertheless, the 
gas is expected to undergo a phase transition to a super- 
fluid state at sufficiently low temperatures. In our opin- 
ion, formulating an accurate microscopic theory [32] for 
this so-called Kosterlitz-Thouless transition is an impor- 
tant challenge for the future. Finally, we would also like 
to mention the possibility of studying mixtures of 
Fermionic and Bosonic atomic gases. One exciting pos- 
sibility is a mixture of doubly spin-polarized atomic ^Li 



and ^Li. In this manner one would create the dilute 
analogue of ^He and ''He mixtures. Such a mixture of 
doubly spin-polarized atomic ''Li and ^Li is rather inter- 
esting because, as we have seen in Sec. 3, the ''Li atoms 
essentially do not interact with each other but only with 
the ^Li atoms. However, due to the latter interaction they 
can now also be evaporatively cooled on a time scale 
which is much smaller than their lifetime. It therefore 
seems worthwhile to extend the methods discussed in 
this paper to this system, and in particular to study 
Bose-Einstein condensation in the presence of an essen- 
tially ideal Fermi gas. 

It is clear from these examples, and from many others 
that we unfortunately are not able to discuss here, that 
we are heading towards an exciting future, in which no 
doubt many surprises will await us. We hope that this 
paper will be of some use to people that intend to con- 
tribute to this rapidly developing area of atomic physics. 
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